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£>x in the chiral limit within the 1/N C expansion. * 
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I report on a recent calculation done in collaboration with E. de Rafael @ of the invariant B K factor of K°-K° 
mixing in the chiral limit and to next— to- leading order in the 1/N C expansion. This calculation is, to the best of 
our knowledge, the first example of a calculation of Bk in which there is an explicit analytic cancellation of the 
/i renormalization scale and the scheme dependence between the Wilson coefficient and the corresponding kaon 
matrix element. I try to emphasize the ideas involved in the approach and how the method could be applied to 
other physical situations, rather than the details of the numerical analysis for which I refer the reader to ref. |?J . 



1. INTRODUCTION 

Real understanding of weak interactions of kaons 
is always hindered by the fact that their QCD in- 
teractions are not perturbative in a s . This is a 
problem with (at least) two well separated scales 
(namely M\y and Aqcd ~ 1 GeV) and, there- 
fore, the well known technique of Effective La- 
grangians can be applied. There is no problem 
as long as one is in the perturbative regime of 
the expansion in a Sl but kaons are lighter than 
Aqcd an d therefore this perturbative regime is 
long gone by the time one gets down to the kaon 
scale. In practice one assumes that the effective 
Lagrangian obtained after integrating the heavy 
fields of the Standard Model (i.e. heavier than 
Aqcd) is valid even immediately below the charm 
threshold. The problem then is to match this 
effective Lagrangian valid at (i < m c to the ef- 
fective Lagragian valid at ^ ~ Mjc, wich is the 
relevant one. 

To be more specific, integrating heavy fields in 
the Standard Model results in an effective AS = 2 
Lagrangian that reads 
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"UAB-FT-496. Talk delivered at the Int. Euro-conference 
on Quantum Chromodynamics, Montpellier, France, Jul 
2000. 



where A„ 



v* d v qs . 



q = u, c, t , with -Fi,2,3 being 
functions of the heavy masses of the fields t, Z°, 
W ± , b, and c which have been integrated out f] 
and 



;as=2 



(as) = [s L (x)^d L {x)\ 



(2) 



with qL = 1 ~ 2 75 q is a four-quark operator. This 
four-quark operator has a hidden dependence 
on the renormalization scale because its ma- 
trix elements are divergent. This divergence is 
subtracted away using the MS renormalization 
scheme, but a fi scale dependence is left behind 
in the process. In fact, because the short distance 
behavior of the operator (J2J) is not the same as the 
matrix elements of the full Standard Model where 
W, etc., are virtually propagating, one must in- 
troduce some coefficients, the so-called Wilson co- 
efficients Cas=2(a0 that correct for this effect. It 
is a consistency condition then that the /i depen- 
dence should drop out of any physical observable, 
since after all one could always do the calculation 
in the full Standard Model (at least in principle) , 
where no /i scale is necessary. Like with the [i pa- 
rameter, the same "declaration of independence" 
also applies to any convention: naive dimensional 
regularization vs. 't Hooft-Veltman, definition of 
evanescent operators, etc... However the ques- 
tion immediately arises, how can kaon matrix el- 
ements of the operator (0) know what you decided 



2 For a detailed discussion, see e.g. refs. 
therein. 
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to do with the commutator of the 75 matrix in the 
calculation of Cas=2 (m) ? For one thing kaon ma- 
trix elements of the operator (||) will depend on 
masses and effective couplings but there can be 
no single Dirac matrix since kaons are bosons, of 
course! The dependence on conventions has to 
be implicitly hidden in what you call couplings 
and massess of mesons when you go away from 
d=4 dimensions in dimensional regularization. I 
shall explicitly demonstrate this point using the 
large- iV c expansion later on. I shall also show the 
explict u cancellation between Wilson coefficients 
and matrix elements of mesons. 

At the scale u ~ Mr chiral symmetry con- 
straints the form of all the possible operators 
appearing in the correspondimng effective La- 
grangian in a chiral expansion in powers of deriva- 
tives and masses. To lowest order, which is 0(p 2 ), 
there is a unique operator with the same symme- 
try properties as Qas=2(x), which is HQ 

9as= 2 (») tr [X 32 (D^)UX 32 (D^)U] ,(3) 

where A32 and U are matrices in flavor space. A32 
denotes the matrix A32 = <5i3<%, and U is the 
3x3 unitary matrix which collects the Goldstone 
fields and which under chiral rotations transforms 
like U — ► VrUVI- The parameter 5as=2(m) is 
a dimensionless coupling constant which depends 
on the underlying dynamics of spontaneous chiral 
symmetry breaking (SxSB) in QCD and that has 
to be determined, as always, by a matching con- 
dition. Leaving aside the known functions i*i,2,3 
it has become conventional to define the so-called 
invariant Br parameter by means of the following 
matrix element 

< K°\Cas=2(^Qas=2(Q)\K° >= p 2 K M 2 K B K .(4) 

Because of the presence of Cas=2{u), Br so de- 
fined must be u (and convention) independent. 
There is also an analogous equation without the 
wilson coefficient Cas= 2 (^) that defines the as- 
sociated Br(u) parameter which, unlike Br, is 
obviously a (and convention) dependent. From 
Eqs. (Q,||) it follows that (with inclusion of the 
chiral corrections in the factorized contribution) 



in other words, Br is also fixed by a matching 
condition. The rest of this paper is a report on 
the calculation of Br carried out in ref . [Q . 

The Wilson coefficient is known from renormal- 
ization group-improved perturbation theory and 
reads 
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Br = -C A s=2(m) x g A s=2(^) , 



(5) 



where the explicit expressions for /3± , P2 , 7i and 
72 can be found in ref. Q and n encodes the 
choices for 75 in d dimensions: k = in naive 
dimensional regularization and n = —4 in the 't 
Hooft-Veltman scheme. Eq. (Q) is obtained af- 
ter insisting on the large- N c limit, which consists 
in taking N c — > 00 with the condition N c a s — ► 
constant J5|. In the strict large- N c limit one finds 

Br = \ , (N c - 00) . (8) 

As to the large- N c expansion, reference || ex- 
plains why taking this limit is a useful thing to 
do because it simplifies QCD without crippling 
it; although the assumption that QCD confines 
in this limit must be made. 

For us it will be very important that, because 
the large--/V c expansion is a systematic approx- 
imation to QCD (both in the perturbative and 
nonperturbative regimes), it can be used as an 
interpolating function between low energies and 
high energies in a systematic way. 

Our calculation will be based on the following 
inputs: 

• The large A^ c expansion. 

• Chiral perturbation theory controls the low- 
energy end of QCD Green's functions. 

• The Operator Product Expansion controls 
the high-energy end of QCD Green's func- 
tions starting from Q ~ lGeV onwards. 

• The curve connecting low-energy and high- 
energy is smooth, i.e. there are no "wig- 
gles" in the intermediate region. Such wig- 
gles would destroy the cherished belief that 
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Nature seems to be "subtle but not mali- 
cious" .0 

Let us now go on with the matching condition 
that 3as=2(a*) must satisfy. From the fact that 
the covariant derivatives contain external fields 
(D^WU = d^WU + iWr»U -il*), one sees that 
Eq. (||) yields a mass term for the external 
fieldQ, among other terms. This r^ s is of course 
the same external field that couples to the right- 
handed current (Ir^sr appearing in the kinetic 
term for the quark in QCD, i.e. 



l- 75 1 + 7e 



q.(9) 



Therefore at the level of the QCD Lagrangian a 
mass term for the r y field can only come about 

as /Q 

because of the operator Qas=2 of Eq. fl2j) in the 
Lagrangian Equating the mass term r r rt 8 
obtained from Eq. (H) to that obtained from Eq. 
(0) (supplemented with £qcd of course) one ob- 
tains the matching condition for g\s=2 ip) R] : 

5AS=2(A*,e)= 1- 

Mhad. 



had. 



y/2 



32ir 2 F 2 



r(2-e/2) 



where /i^ad ^ s a scale used for normalization, e.g. 
z = Q 2 /Mhad > but otherwise totally arbitrary. 
Therefore the apparent quadratic dependence of 
Eq. (|l0|) on /ihad. is actually fictitious. There will 
be nevertheless a certain /ihad. dependence in the 
final answer for Bk due to the fact that there is 
always a residual p dependence because the effec- 
tive Lagrangian of Eq. ([!]) is constructed only to 
next-to-leading order in a s (see Eqs.(^,^)) and we 
will identify p = /ihad. as a typical hadronic scale. 
For the final numerical applications we shall take 
1 GeV < ^had. < 3 GeV. 

The divergent (i.e. ill-defined) integral defin- 
ing 5as=2(Mi £ ) m Eq. ( [l0| ) must be under- 
stood within the same regularization procedure 
and conventions used to compute Cas=2(^) in 
Eq. (0), i.e. MS and the definition of evanescent 
operators of ref. || . 

3 There's no free lunch. Any interpolating procedure needs 
two good ends to bridge from, and the assumption that 
nothing wild is going on in between !! 
4 Flavor indices are important ! 



The function W(z) is defined as 
W(z) = ^<] lfi ( Z i.) , 



(11) 



where Wlrlr ^ s defined through an integral over 
the solid angle of the momentum q, 



unfactorized 
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and Q 2 = —q 2 . In Eq. ( |T^ ) the function 
W lrlr{i, stands for 
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(13) 



(14) 



(15) 



Since W(z) ~ C(N°) and F 2 - 0{N C ) one sees 
that Eq. @ yields ffAS=a&*,e) = 1 - 0(1/N C ) 
with the subleading term governed by W(z). 

What do we know about W(z) ? In general 
W(z) is a complicated function of the variable z. 
However in the large- ./V c limit it simplifies some- 
what and can be expressed as an infinite sum of 
single, double and triple polesQ in the form 



W[z\ = 6 



£7 



■E 

1=1 



p 2 
i=l Hl 



(z + pi) (z + pi) 2 {z + pif 



,(16) 



where pi = M 2 //ihad an( l ^ s ^ ne mas s of the 



«-th zero- width resonance. If one knew all the 
masses and constants appearing in W(z) from the 
solution to large-iV c QCD one could use this in- 
formation, plug it in Eq. (10) and compute Bk 



5 J. Bijnens and J. Prades are thanked for pointing out 
the existence of triple poles, which were actually ignored 
in the previous version of this analysis. See Ref. M. 
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through Eq. (||). Since this solution is not avail- 
able, we will have to follow a different strategy. 
Firstly, one recalls that it is due to the pres- 
ence of a perturbative contribution to a given 
Green's function that one must take these sums 
over resonances truly up to infinity because oth- 
erwise parton-model logarithms cannot be repro- 
duced j6|. Clearly, when the sums extend all the 
way to infinity, in general one cannot expand in 
powers of 1/z at large 2 in a naive way. However 
our function W(z) is an order parameter of spon- 
taneous chiral symmetry breaking or, in other 
words, the contribution of perturbation theory to 
W(z) vanishes to all orders. One can then try 
to match a naive expansion in 1/z to the Opera- 
tor Product Expansion considering finite sums as 
an approximation to the originally infinite sums. 
This way one can fix the unknown parameters and 
masses through the knowledge on the behavior of 
W(z) at large z given by this Operator Product 
Expansion. 

We have no proof that these infinite sums have 
to behave as if they were finite but we know 
that, if this is so, a successful framework such 
as Vector Meson Dominance then becomes the 
leading term in the approximation to large N c 
that wc propose, instead of a phenomenological 
ansatzQ. Besides, our approximation works well 
in those cases where we can compare with the 
experiment [ ^0|Jfl| ] . 

Consequently we impose the following behavior 
on W(z): 

• i)At low z, from chiral perturbation theory, 
one finds IT^ITT _ 



W(z) 



6 







5L 2 + L 3 + L 9 ) z (17) 



Notice that in i) chiral logarithms are neglected, 
in accord with the large--/V c expansion, ft is a 
phenomenological fact that they are small in vec- 
tor and axial-vector channels, which are precisely 
the relevant ones here. For the Li couplings we 
use lowest meson dominance, as explained in ref. 
fl(],[|- This also works phenomenologically. The 
factor of 6 in Eq. ( |l7| ) explains the corresponding 
factor in W{z) in Eq. |L6|). 

As to ii) , notice the e and k dependence of the 
result. This is brought about by the need to do 
the calculation in d dimensions, as repeatedly em- 
phasized. Then the appearance of structures like 
7a(1 -75)7p7a*®7 A (1 -7s)7 p 7 m in the OPE cal- 
culation of W(z) makes the definition of 75 and 
evanescent operators of practical relevance. We 
followed closely refs. 

Considering the higher orders in Eq. ( p"8| ) there 
has recently been the warning, issued by the au- 
thors of ref. Jig ], that these subleading effects 
in the OPE expansion might actually be numeri- 
cally important. 1 don't expect this to be the case 
here because, as we shall see, the impact of the 
OPE fall-off in the final numerical value is not 
large, but a more detailed analysis is currently 
underway fl6||. 

All in all, Eqs.(|l7],[lJ) impose the following con- 
ditions on the parameters and masses of W(z): 
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from 0(z) terms in Eq. (p7[). Furthermore, 
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from the absence of 0(1/ z°) terms in Eq. 
also 



(19) 
(20) 



ii) At high z, from the operator product ex- 
pansion, one getsQ 



W(z) « 
24^a s F 2 

Mhad. 
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(21) 



(22) 
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Figure 1. Plot of the hadronic function W[z] in Eq. versus z — Q 2 /^ ad for ^had. = 1-4 GeV, 

Fq = 85.3 MeV and My — 770 MeV. The red curve is the function W[z] corresponding to the contribution 
of a single p resonance (see text). The green line represents the low-energy chiral behaviour and the blue 
curve the prediction from the OPE. 
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from the 0(l/z) term in Eq. 

When all these constraints are imposed on 
W(z) in Eq. ( |l6| ) the expression for <?as=2(m) 
in Eq. ( jlpj ) in the MS renormalization scheme 
reads 



5as=2(m) = 1 - 



Phud. 

32tt 2 F 2 



ft log 



Phad. 
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C*s(Mhad.) 

OC 
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-a, log pi H h - — 

Pi 2 pi 



(23) 



where in the second line we have written 
the renormalization group- improved result. We 
would like to insist on the fact that this result, 
contrary to the various large-iV c inspired calcu- 
lations which have been published so far 17 1^, 
is the full result to next-to-leading 
order in the 1 /N c expansion. The renormalization 
p scale dependence as well as the scheme depen- 
dence in the factor S cancel exactly, at the next- 
to-leading log approximation, with the short dis- 
tance p and scheme dependences in the Wilson 
coefficient Cas=2{p) m Eq. (Q). Therefore, the 
full expression of the invariant -Bi<— factor to low- 
est order in the chiral expansion and to next -to- 
leading order in the l/N c expansion, which takes 
into account the hadronic contribution from light 
quarks below a mass scale /ihad. is then 
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a s (^had.) 
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Mhad. 
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111 
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.(24) 



As remarked earlier, the numerical choice of 
the hadronic scale /ihad. is, a priori, arbitrary. 
In practice, however, /ihad. has to be sufficiently 
large so as to make meaningful the truncated 
pQCD series in the first line. Our final error 



in the numerical evaluation of Bk will include 
the small effect of fixing /ihad. within the range 
1 GeV < /i had . < 3 GeV. 

How many resonances should one include in the 
sum in Eq. (24)? Clearly, the more the bet- 
ter. However, given that there is only a limited 
amount of information coming from Eqs. (17, 18]) 



one can in practice determine the couplings of just 
one resonance, the p vector meson. With them, 
and given the rho mass, the corresponding W(z) 
function of Eq. ( |l6| ) is shown in Fig. 1, where 
it is seen that W(z) intersects the OPE curve at 
a certain value of z that we will call z. It turns 
out that z ~ 0.39 (i.e. Q ~ 900 MeV). As can 
be seen, the matching to the OPE curve is not 
perfect and this is why, for the final numerical 
estimate, we propose to use the hadronic curve 
W(z) only in the region < z < z while in the 
range z, < z < oo we use the first term of the 
OPE, Eq. (|l8|). This results in a different expres- 
sion for the Bk parameter which reads 
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Of course in the limit z — * oo, Eq. ( |25| ) repro- 
duces Eq. (^J) , but this limit in practice can only 
be achieved by using more and more resonances 
in the sum, which requires the knowledge of more 
and more terms in both the low and high-energy 
expansions of W(z). The physical interpretation 
of the result ( p^ ) is clear: the area under the curve 
in Fig. 1 is roughly the amount one subtracts 
from the large- N c limit result of Bk = 3/4. In 
numbers one gets 



B K = 0.38 ±0.11, 



(26) 



where the error is the combined limit of errors 
from short- (e.g. how much is and long- 

distance (e.g. how much is the p mass?) contri- 
butions to Bk for a choice of /ihad. in the range 
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1 GeV < /ihad. _ 3 GeV. I refer the reader to ref. 
for details. This result, of course, does not 
include the error due to next-to-next-to leading 
terms in the 1/N C expansion, nor the error due 
to chiral corrections in the unfactorized contri- 
bution, but we consider it to be a rather robust 
prediction of Bk in the chiral limit and at the 
next -to- leading order in the 1/N C expansion for 
the following reasons. Firstly, Fig. 1 shows that 
most of the area under the curve is in the region 
< z < 0.4 which makes the final number for 

Bk relatively insensitive to the part of the curve 
coming from the OPE. Secondly, at a more qual- 
itative level, if one accepts that the low-z region 
is correctly given by Chiral Perturbation Theory 
up to a scale of z ~ 0.06 (i.e. E ~ 350 MeV) 
and the OPE also describes reasonably well the 
high-z region starting from z ~ 0.5 (E ~ 1 GeV); 
then the conclusion is that there isn't much room 
left in between to do a lot of things ! 

Our final result in Eq. ( |2^) is compatible with 
the current algebra prediction |20| obtained from 
the K + — > 7r + 7r decay rate at lowest order in chi- 
ral perturbation theory. It turns out that the 
bosonization of the four-quark operator Qas=2 
and the bosonization of the operator Q2 — Qi 
which generates AI =1/2 transitions are related 
to each other in the combined chiral limit and 
next -to-leading order 1/N C expansion. It fol- 
lows that decreasing the value of Bk from the 
large- N c prediction of 3/4 down to the result in 
Eq. (|2^) is correlated with an increase of the cou- 
pling constant in the lowest order Effective Chi- 
ral Lagrangian which generates AI =1/2 tran- 
sitions, providing therefore a first glimpse at a 
quantitative understanding of the dynamical ori- 
gin of the AI — 1/2 rule. In the future we would 
like to make some refinements and a detailed com- 
parison between our result and other previous es- 
timates [jl6). Among these "refinements" one has 
most notably the inclusion of higher order terms 
in the chiral expansion. These are crucial in order 
to be able to make a meaningful comparison with 
the values favoured by lattice QCD determina- 
tions |22| as well as by recent phenomenological 
analysis p3| , p4[ . 

I would not like to finish without stressing 



that our framework can in principle be improved 
systematically, for instance by computing more 
terms in Eqs. (0,|l|). It can be applied as well 
to other physical situations involving weak ma- 
trix elements. In order to do this one might find 
helpful the following 

DO-IT-YOURSELF KIT: 

• i)Pick the QCD Green's function, Q, which 
defines the matrix element you are inter- 
ested in. 

• ii)Find the large-_/V c resonance representa- 
tion of Q (i.e. the interpolating function). 

• iii)Impose on the above interpolating func- 
tion the long- and short-distance con- 
straints stemming from Chiral Perturbation 
Theory and the Operator Product Expan- 
sion, respectively. 

• iv)Compute number. 

• v)Enjoy. 
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DISCUSSION 

H. Fritzsch (Munich) 

I would expect that you encountered problems 
with your technique if the singlet h is involved 
since the role of the gluon anomaly is different for 
large N c . 

S.P. 

As is well known, the gluon anomaly is at the 
origin of the rj mass. In the strict large N c limit 
the anomaly is switched off and M n ' — > 0. It may 
seem then that the large- N c expansion is a very 
bad expansion for the 7/ since in the real world 
M v i ~ Mproton ~ 1 GeV >> M n ; i.e. it doesn't 
look at all like a Goldstone boson which is what 
the N c — 00 result says. 

Although I don't know whether the 1/N C ex- 
pansion will be in the end really helpful for cal- 
culations involving the h channel, I would like 
to explain why I think that the above argument 
is potentially fallacious. Take good old QED as 
an expansion in powers of a around a — 0. If 
you were to compare the real world (where e - 
and e + certainly scatter) to what you find from 
the strictly speaking first term in this expansion 
, which is the one at a — (i.e. no scattering), 
you could also conclude that the a expansion is 
a disaster. Well, in some respects a — corre- 
sponds to N c = 00. In other words, one should 



always go to first nontrivial order in the expan- 
sion before comparing to the real world. It is not 
unconceivable then that, after the bulk of the rf 
mass is included in the first nontrivial term in 
the 1/N C expansion, the rest of the series in 1/N C 
turns out to be a well-behaved innocent-looking 
power series. 

Having said this, I also want to stress that in 
our particular case of Bk, at the order we are 
computing, the h channel doesn't play any role 
at all. It's all governed by vectors and axial- 
vectors for which we know that resonance sat- 
uration (which is basically what the large- N c ex- 
pansion amounts to in this case) is a good thing 
to do. (I thank M. Knecht for reminding me of 
this.) 



